LAMMPS (10 Jan 2012)
# LAMMPS benchmark of granular flow
# chute flow of 32000 atoms with frozen base at 26 degrees

variable	x index 1
variable	y index 1

units		lj
atom_style	sphere
boundary	p p fs
newton		off
communicate	single vel yes

read_data	data.chute
  orthogonal box = (0 0 0) to (40 20 37.2886)
  2 by 1 by 2 MPI processor grid
  32000 atoms
  32000 velocities

replicate	$x $y 1
replicate	2 $y 1
replicate	2 2 1
  orthogonal box = (0 0 0) to (80 40 37.2887)
  4 by 1 by 1 MPI processor grid
  128000 atoms

pair_style	gran/hooke/history 200000.0 NULL 50.0 NULL 0.5 0
pair_coeff	* *

neighbor	0.1 bin
neigh_modify	every 1 delay 0

timestep	0.0001

group		bottom type 2
3648 atoms in group bottom
group		active subtract all bottom
124352 atoms in group active
neigh_modify	exclude group bottom bottom

fix		1 all gravity 1.0 chute 26.0
fix		2 bottom freeze
fix		3 active nve/sphere

compute		1 all erotate/sphere
thermo_style	custom step atoms ke c_1 vol
thermo_modify	norm no
thermo		100

run		100
Memory usage per processor = 37.0464 Mbytes
Step Atoms KinEng 1 Volume 
       0   128000    3136556.5    6404.5051    119323.52 
     100   128000      3137160    6284.0549    119327.22 
Loop time of 0.893496 on 4 procs for 100 steps with 128000 atoms

Pair  time (%) = 0.524706 (58.7251)
Neigh time (%) = 0.0462004 (5.17075)
Comm  time (%) = 0.0590464 (6.60847)
Outpt time (%) = 0.00043273 (0.0484311)
Other time (%) = 0.26311 (29.4473)

Nlocal:    32000 ave 32000 max 32000 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Nghost:    5456 ave 5456 max 5456 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Neighs:    115135 ave 115135 max 115135 min
Histogram: 4 0 0 0 0 0 0 0 0 0

Total # of neighbors = 460540
Ave neighs/atom = 3.59797
Neighbor list builds = 2
Dangerous builds = 0
